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Abstract 

By using the extension of the statistical Wick's theorem (Gaudin's theorem) to deal with gener- 
alized statistical density operators (those which can be expressed as the product of and operator 
carrying out a canonical transformations times a density operator) and using the appropriate limits 
we are able to rederive in a very simple way the standard generalized Wick's theorem for overlaps 
of mean field wave functions. Due to the simplicity of the derivation it is now straightforward to 
consider more involved cases and some of them are discussed. The present derivation also allows 
to obtain general and compact formulas for other particular cases of the generalized Wick theorem 
involving overlaps of multiquasiparticle excitations of product wave functions. The new expressions 
allow to reduce the combinatorial complexity of the standard calculation of the above overlaps. 
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I. INTRODUCTION 



Nowadays it is widely recognized the need for going beyond the mean field approximation 
in order to better understand the low lying structure of the atomic nucleus (see Ref [1] for a 
thorough discussion of the topic). The mesoscopic nature of the nucleus makes essential the 
restoration of the symmetries broken at the mean field level in order to incorporate into the 
wave function the expected behavior under the symmetry's operators that are responsible 
for the selection rules of many observable quantities. The correlation energy gained by the 
restoration of symmetries and the consideration of quantum fiuctuations in relevant degrees 
of freedom plays a very relevant role not only in the description of the nucleus' binding 
energy but also in the distribution of states in the spectrum of the nucleus. In the process of 
incorporating into the theoretical framework the quantum fiuctuations mentioned previously, 
tools like the generalized Wick's theorem are widely used. The generalized Wick's theorem 
(GWT) for fermions allows the evaluation of the overlap of a general product of fermionic 
creation and annihilation operators between any two non orthogonal product wave functions 
|$o) and According to the GWT, the overlap ($i|0|$o)/(^i|^o), where O is a general 
product of creation and annihilation (5k quasiparticle operators, is given as the sum of 
all possible contractions of all possible combinations of two quasiparticle operators, namely 

Cii mn = ($l|/3r„/3n|^o)/($l|$o) 

Ci2 mn = mPml3i\^,)/m^0) (1) 

C21 mn = ($l|/?+/?„|$o)/(*l|*o) 

C22 mn = ($l|/3+/3+|$o)/($l|$o) 

These four contractions can be gathered together as members of the bipartite matrix of 
contractions 

cJ^J^^A (2) 

\ C2I C22 / 

that will prove to be a convenient quantity in future developments. Both the proof of the 
generalized Wick's theorem as well as the explicit expression of the contractions have been 
given by several authors [2-4] using different approaches. In all the cases the proof of the 
theorem as well as the contractions are obtained after quite lengthy considerations and in a 
manner that is not prone to generalizations. On the other hand, the generalization of the 
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standard Wick's theorem for mean values to the framework of an statistical admixture of 
quantum states (where mean values are replaced by traces over the Fock space including 
a statistical density operator) is rather straightforward and the proof of the theorem is 
fairly simple [5] owing to the simplifications introduced by the cyclic property of the trace. 
The generalization of this theorem (Gaudin's theorem thereafter) to the case where the 
statistical density operator is replaced by a more general operator given by the product of 
the density operator times an operator carrying out a canonical transformation is also rather 
straightforward [6, 7]. The reason is that the spirit of the proof of Gaudin's theorem (the 
cyclic invariance of the trace) is still valid and the proof of the generalized theorem proceeds 
along the lines of the standard one. The purpose of this paper is to show that the generalized 
Wick's theorem can be deduced from the corresponding statistical version by an appropriate 
limiting procedure of the "probabilities" of the statistical admixture. This limiting procedure 
together with the simplicity of the statistical version of the theorem will allow us to derive 
formulas for the most general situation in a fairly simple way. By generalizing the previous 
limiting procedure we will be able to compute not only the standard overlaps but also 
overlaps of multiquasiparticle excitations that can be rather cumbersome to compute by 
using the standard techniques [8, 9]. The novel expressions obtained in this paper for the 
contractions are very compact and therefore their use substantially simplifies the usually 
lengthy calculations (with a complexity of combinatorial type) involved in considering those 
multiquasiparticle overlaps. 

II. DEFINITIONS AND PROOF OF GAUDIN'S THEOREM 
A. Basic definitions 

First of all let us introduce some basic notation that will be used to simplify the expres- 
sions derived below (the notation closely follows that of [2]). A given set of fermionic creation 
and annihilation quasiparticle operators (3k and (3^ defined in terms of a reference single 
particle creation and annihilation set of operators bi and bi by means of linear combinations 

Pk^j:iA;,bi+B:,bi 

= J2i ^ikbi + Dikbi 



3 



The reference single particle creation bi and annihilation bi set of operators is assumed to 
satisfy canonical anticommutation rules {bi,bii} — {bi,bi/} — and {bi,bii} — but the 
operators bi and bi (as well as Pi and Pi) are not necessarily related by hermitian conjuga- 
tion. In the following, creation and annihilation operators expressed with Latin letters will 
correspond to single particle operators whereas the ones in Greek letters will correspond to 
quasiparticles. For both kinds of operators Latin subindices will be used. These and other 
relationships to be considered below can be written in a more compact way by introduc- 
ing vectors of dimension 2N [N is the total number of creation or annihilation operators 
considered) incorporating both annihilation and creation operators both for the quasiparticle 

= (Pi, P2,---, Pn, Pi, P2,---, Pn) (4) 

and the single particle 

= {bi,b2, . . . ,bN,bi,b2, . . . ,bN) (5) 

operators. For these vectors of operators Greek indices like the n above are introduced. It 
is also convenient to introduce the notation 

(Pl,p2,---,PN,PI,P2,---,PN) (6) 

as well as the 2N x 2A'" matrix 

/ n 1 \ 

(7) 




that allows us to express the ct^u set in terms of and viceversa as — ^j^cr^^a^. Now 
the linear combination of Eq. (3) can be written in a more compact way as 



a, 



1/ 

where the 2N x 2N matrix W has the following bipartite structure 

(9) 




If it is required the new set of quasiparticle operators to satisfy the canonical fermionic 
commutation relations 

{a^, a^} = 6 1^1, (10) 
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(and assuming that the single particle set does) then the coefficients W of the linear 
combination defining in terms of must satisfy 



which is the condition for the transformation matrix to be a canonical transformation. 

Nothing has been said about hermiticity since the operators bk need not to be the hermi- 
tian conjugate oibk (the only relevant property are the anticommutation relations). However, 
and in order to simplify the following considerations we will assume in the following that 
the single particle operators 6^ = bk are the hermitian conjugate of the bk. Concerning 
the quasiparticle operators (3i we will assume that they are not in general the hermitian 
conjugates of (3i. However, if we require (3i = to be the hermitian conjugate of I3i then 
the coefficients of the transformation relating them to the set bk and b'l must satisfy the 
two requirements A — D and B — C. In terms of the W coefficients this translates to 
aWa = W*. 

In order to prove the statistical Wick's theorem (or Gaudin's theorem) mentioned in 
the introduction we only have to assume that there is a statistical density operator p, not 
necessarily hermitian, and satisfying 



where the matrix M is in principle arbitrary. However, if this transformation is applied to 
the canonical commutation relation p^^ {a ^, a^} p = a^y we obtain {p^^ct^p, p^^ctj,/)} = a^^, 
and using Eq. (12) the condition McrM^ = cr is obtained for the matrix M implying that this 
matrix indeed represents the one of canonical transformation (see Eq. (11)). Usually, the 
density operator as well as other operators carrying out canonical transformations and to be 
used later are written in terms of the exponential of one body operators (see Appendix A). In 
this case the transformation matrix M is given by the following exponential M — exp(— crK^) 
where is a skew-symmetric matrix (see appendix A for details). However, not all the 
matrices satisfying the canonical transformation condition of Eq. (11) can be written as the 
above exponential (the most general transformation in a space of dimension 2 is enough to 
find a counterexample, see appendix A) and in most of the exposition we will consider the 
density operator p as a general operator not necessarily expressible as the exponential of an 
one body operator. 



(11) 




(12) 



V 
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We also have to introduce the concept of the trace of an operator A over the whole Fock 
space 

Tr[i] = {<P\A\<P) + Y,WmA(3M + ^ J](0|/3^/3ni/5+/3+|0) + • • • (13) 

m mn 

where |0) is a product wave function and /3m and (3^ are the corresponding annihi- 
lation and creation operators associated to The trace satisfies the cyclic property 
Tr[AiA2 ■ ■ ■ Ak-iAk] = Tr[AKAiA2 ■ ■ ■ Ak-i] that will turn out to be the fundamental 
property to prove Gaudin's theorem. 

B. Proof of Gaudin's theorem 

Now let us consider the evaluation of the following trace 

] (14) 

with K an even integer number. The result to be obtained is nothing but Gaudin's theorem 
of Ref [5] but it is repeated here as it will be instructive for further generalizations to 
be considered below. We start by moving to the right by using the anticommutation 
relations 

Tr[pa^^an2 ■ ■ ■ a^^] = {a^j, a^jlTrfpa^gQ;^^ • • • a^^] - Trfpa^^a^^ • • • a^^] (15) 

and keep moving it to the right until we reach the last quasiparticle operator in the argument 
of the trace 

Tr[pQ!jitjQ!^2 • • • Q;^^] — {q;^]^, (y,^2Y^\p^iJ,z^iJ,A ' ' ' ^ij-k\ 

{ojjui, Q^//3}Ti'[pQ;/i2'-'^/i4 ' ' ' '-'^/i/f] 
+ h 

- Tr[pQ;^2Q;^3 • • • a^^a^ij (16) 

Using now the cyclic invariance of the trace the last term is written as Txlan^pa^^an^ ■ ■ ■ a^^] 
and using the property of Eq. (12) we arrive at 

] = ^M^,,,Tr[pa ] (17) 
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Moving now this term to the left hand side of Eq. (16) we obtain 

+ M^^^jTr[pQ;^^Q;^2 • • • a^^] ^ {a„^, ai^^}Ti[pai^^a^^ ■ ■ ■ a^^] (18) 

Ml 

- {tti.!, a^3}Tr[pa^2"M4 " " • "mx] 

that can be written as 

Tr[pQ;^jQ;^2 • • • a^^] = C^^^jTrfpa^gQ;^^ • • • a^^] - C^^^gTrfpa^^a^^ • • • a^^] H (19) 

where the matrix of contractions C^^^j have been introduced. They are given by 

C^,,, = 5](1 + M);l^{a,^, a,,} = [(1 + M)-V]^,^, (20) 

Now repeating the above procedure for the remaining traces of Eq. (19) we arrive at the 
Gaudin's theorem which states that the trace 

IV[p] ^^'> 

equals the sum of the product of all possible contractions of Eq. (20). It is important to 
point out here that the above derivation is independent of the product wave function and 
quasiparticle operators entering the definition of the trace in Eq. (13) and the only relevant 
properties are the cyclic invariance of the trace and the transformation law of Eq. (12). 



C. Operator overlaps as limits of statistical traces 

Usually, the statistical density operator p is given as the exponential p — exp{—K) of an 
one-body hermitian operator K — + | ^^^^ K^i^a+ai, with 



(22) 



The matrix K^^ is in this case hermitian whereas K"^^ is skew-symmetric. The matrix M 
corresponding to the representation of the statistical density operator in the quasiparticle 
basis is given by M = exp(— K) and by working in the basis where K is diagonal 





(23) 



the matrix M becomes diagonal 

^ Pi 



(24) 



I/Pi 

with pi = exp{—ki). We will introduce the quantity — 

(Pi)P2, • • • ,PAr, 1/pi, l/p2, • • • , I/pat) such that M^,^ = p,j,5^y. The statistical density 
operator can be now written as 

P = + Y.PmPMWm + ^ 5^PmPn/3+/3+|0)(0|/3n/3^ + • • • (25) 

m mn 

where the quasiparticle operators and are the ones for which the operator K (and the 
matrix representation K) are diagonal and the product wave function |0) is the corresponding 
vacuum of the annihilation operators We also notice in the previous expression of the 
statistical density operator that the quantities Pm can be interpreted as the "probability" of 
a given quasiparticle excitation. By using the previous expression of the statistical density 
operator we obtain 

Tr[pi] = {<i>\m + Y,Pm{<l>\PmAf3M + I J]p„.Pn(0|/5„^/5ni/5+/3+|0) + ■ ■ ■ . (26) 

m ' mn 

that will prove to be useful in the following. 

Let us now consider the generalization of Gaudin's theorem to the case where the statisti- 
cal density operator p is replaced by the product pT~^ where T is an operator performing a 
canonical transformation in the quasiparticle operators cc^ such that the transformed quasi- 
particle operators 

a^^fa^f-' ^J2T^,a, (27) 

HI/ 

still satisfy the canonical commutation relations but not necessarily hermiticity. The re- 
quirement of preserving the canonical commutation relations implies TaT'^ — a. We now 
have 

Y^pf-' = J2 P^T^,pf-'a, (28) 



a,. 



or M^j, = p^T^i, and we are in the general situation that allowed us to derive Gaudin's 
theorem. Therefore, for the evaluation of the trace 

(29) 
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we can use Gaudin's theorem with the contractions given by Eq. (20) and with M = p^T^i^ 
replacing the matrix M. Additionally, we have that 

l^[pf-'A] = {^\A\(t>) + + ^ J]p„p„(^|^^^„i/?+/?+|</)) + • • • . (30) 

m ' mn 

where (0| — (0|T~^ and — 'T'PmX~^- In the limit — > we have 

*M = Mm ly^. (31) 

(010) Tr[pr-i] 
This expression is the link between the results obtained with the statistical density operator 
and the ones obtained with pure states and it will allow to prove in the next section the 
extended Wick's theorem out of the Gaudin's one. 



III. PROOF OF THE EXTENDED WICK'S THEOREM 

In the previous section we have just shown that the overlap between different mean field 
wave functions of a product of quasiparticle operators can be written as the limit of the 
corresponding trace when the statistical probabilities go to zero 

(010) Tr[pr-i] 
The argument of the limit of the right hand side of this equation can be easily computed 
using Gaudin's theorem. It is given as the sum of all possible contractions between pairs of 
quasiparticle operators that are given by 

C^,^, = (1 + M)-l,^ = [(1 + M)- V]^,^, (33) 

The relation given by Eq. (32) tell us that the overlap {(t>\oi^j^^a^^ ■ ■ ■ q;^^|0)/ (0|0) is given as 
the product of all possible contractions C^^j^j between pair of quasiparticle operators which 
are given by 

C^,^, = lim C^iM2 = [(1 + M)-V]^,^, = lim [(1 + pT)-'aU^, (34) 

Pm— »0 Pm— >0 Pm— >0 

provided that the limit exists. In order to compute the limit we have to consider the general 
bipartite structure of the matrix T 



T = (35) 




in order to write 

(1 + pT)-^ 




Y X I \ p 
This expression allows a trivial evaluation of the limit 



lim (l + pT)-^ 



1 

Y X 





that leads immediately to 

C 




or 



m 



mim2 



m 



mim2 



(36) 



(37) 



(38) 



(39) 
(40) 
(41) 



which is the expected result of the GWT (see Refs. [2-4] for details). Another way to 
perform the limit is to use the property of the matrix T (consequence of being the matrix 
of a canonical transformation) 

TaT^a = 1 (42) 
in order to write = aT'^a. Using this property it is very easy to show that 

(1 + pT)-i = (TT^a((TT^(T + p)-^ (43) 



This expression will prove to be useful in taking the limit as it can be easily shown that 

-1 

p \ I A B 



lim 

Pm^O 



1/p 



+ 



C D 




(44) 
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Using now that 



aT^a= 45 





we easily arrive at 

;im .T-.(.T-. + p)-' = I : I = I ...1 . : I (46) 

which is equivalent to Eq. (37) because as can be deduced from the 

property of Eq. (42). 

As a demonstration of the usefulness of the method just described we will consider the 
generalized Wick's theorem for the quantity 

(0l|02) 

where the quasiparticle operators a^^^ are not related to the mean field wave functions |0i) 
of the previous overlap. The above overlap can be written as the following limit 

where it has been assumed that there exist operators % relating the mean field wave functions 

with the vacuum of the quasiparticle annihilation operators of the generalized set 
(/?fe|0) = 0) through the relation = Considering the extended Gaudin's theorem 

of the previous section we can write the argument of the limit of the right hand side as the 
sum of all possible contractions 

<^,.,. = (1 + VTi);i, = {T-,\T,T-,' + P)-'T,)^_^ = {aTla{T,aTla + pY'T^o) ^^^^ 

(49) 

The limit — >^ in this expression can be evaluated straightforwardly by using the result 
of Eq. (44) and the notation of Eq. (35) 

C ^ C J U ^ (50) 

y Yl U( j \ Q j \X2Y2 j 

Finally, we shall derive all the contractions needed for the evaluation of the most general 
overlap 
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(0i|t3|02) Trfptf^tsta] 

(51) 

In this case, the operator 7^ separates the product of quasiparticle operators in two groups, 
the ones to the left of this operator and the ones to the right. Because of this separation we 
will need to consider three different kinds of contractions depending upon in which group 
are located each of the two quasiparticle operators involved in the contraction. To be more 
specific, we will need the contractions 

C2 = ((1 + T3%-VTi)-V)^^ (52) 

when both indices /J, and are in the set of operators to the left of and with indices in 
the set {/ii, . . . , /ij}, 

eg] = ((1 + T3-^T,-VTi)-iT3<7)^^ (53) 

when IJ, is in the set {pi, . . . , Hj} (i.e. to the left) and i/ is in the set {/Jj+i, • • • , I^k} (i-e. to 
the right) and finally 

CS-((1 + T2VTiT3-')-V)^, (54) 

when both indices /i and u belong to quasiparticle operators that are to the right of % and 
therefore are in the set {fJj+i, • • • , px}- The limits ^ in the above contractions can be 
very easily performed using previous considerations like the ones leading to Eq. (50) and 
will not be given here. 

IV. EVALUATION OF THE PARTITION FUNCTION OR NORM OVERLAP 

In this section we will evaluate the expression for the partition function or norm over- 
lap appearing in the previous section. A typical partition function to evaluate is of the 
form Tr[/57^^^]. Taking into account that the density operator also performs a canonical 
transformation, the product pT^^ also corresponds to a canonical transformation whose 
matrix representation is the product of the matrix representation of p and There- 
fore, we could reduce the calculation of Tr[p7^~"^] to the one of Tr[T] where p7[~^ — T. 
However, in the calculations of the traces we will have to assume that at least one of 
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the canonical transformations can be written as the exponential of an one body opera- 
tor and therefore the evaluation of TrfpTj"^] can be considered as more general by assuming 
that the transformation Tj"^ is given by the exponential of an one body operator, namely 
Tf^ = e^^'^ = cxp(— i '^j^) ^^11*^1 considering that p is a general operator per- 

forming a canonical transformation and not necessarily expressible as the exponential of an 
one body operator (see appendix A). To compute Tr[p7i~^] we consider the function 

n(A) = Tr[pexp(-A^i)] (55) 

that reduces to the quantity we want to compute in the limit A = 1. According to previous 
notation the operator 5*1 is given by 5*1 = | YI^v'^ij. <^i^- To evaluate n{X) we consider 

its derivative with respect to A 

= -Tr[pexp(-A5i)^i] (56) 

that can be easily computed with the help of the extended Gaudin's theorem considered in 
the previous section 

^ = -^E(^iUTr[pexp(-A^O«,a.] = -^^(^0,. ((1 + MT(A))-V)^^ n(A) 

(57) 

where, as usual, the transformation matrix T(A) is given by 

T(A) = exp(-A(j§i,A) (58) 

with = - Sf) and M is given by Eq. (12). Using the property (1 + MT(A))-V + 
(7(1 -|- T^(A)M^)~^ = (7 deduced in appendix A and a little of algebra we obtain 

= + MT(A))-'] - iTVl.S,l (59) 

The first term of the right hand side of the equation can be written (see appendix B) as 

Tr[cr§i,^(l + MT(A))-^] = :^Tr[ln(l + MT(A)] (60) 
which allows to integrate the differential equation defining n{X) to obtain 

n(A) = Ce-^^l^^il[det(l +MT(A))]^/2 (61) 
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where C is an arbitrary constant and the relation det A = exp(Tr[lnA]) has been used. The 
constant C is determined by taking into account that n(0) = Tr[p] = C[det(l + M)]^/^ so 
that we finally arrive to the general expression 



Tr[prr^] ^ [det(l + MT(l))]V^ 

Tr[p] [det(l + M)]V2 



(62) 



In the very common case where p is also given as the exponential of an one body operator 
we can choose p — Ti and deduce that 



Tr[p] =e^^[^^''l[det(l + M)]V2 



(63) 



where we have made use of the property MT(1) = 1 holding in this specific case. Taking 
everything into account we arrive at the particular result 



TrfpTf ^] = e-^^['^(^i-^^)l[det(l + MT(1))]V2. 



(64) 



In the case where p is "diagonal" in the sense that it is given by 

p = exp(^ km&m) = exp(- ^ k^) exp(- ^ (S^)^^ a^) (65) 



HI/ 



with 



-k 
k 



(66) 



the transformation matrix M is diagonal and given in terms of the probabilities Pk and their 
inverses 

p 



exp(cT§p) =1 I = P- 

p-^ 

The constant term in Eq. ( 65) is written in terms of the probabilities as 



(67) 



exp 



1/2 



Pn 



det 



1 

p 



1 1/2 



(68) 



Using all these definitions we can finally write 

1/2 



(69) 
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that leads, by means of Eq. (64), to 







det 


:)1 







1/2 



which is finally written as 

Now the p ^ limit is straightforward 



= hm TrfpTf 1] = g-^^t'^^il 



det I I ^ M + I ^ ° I Ti 
Op/ \ 1 



-I 1/2 



[det(l+pTi)]i/2 
(70) 

(71) 




1/2 



= e-^^['^^il[detXi]^/' (72) 



and corresponds to the expected result [2-4]. 

V. MULTIQUASIPARTICLE OVERLAPS 

Finally, we will consider a very helpful extension of the above considerations to the 
evaluation of overlaps between multiquasiparticle excitations 







(73) 



where the "tilded" quasiparticle operators are related to the "untilded" ones through a canon- 
ical transformation 

V = ra^r-i = ^T^,Q;, (74) 



a. 



and \4>) = T\4>) . In the spirit of the preceding sections the above overlap can be written as 

(75) 



where the limit has to be understood as a limiting process where the probabilities p^^ , • • • , Pjxk 
are made to tend to infinity while at the same time the other probabilities {u ^ fj,i, . . . , ijlk) 
are set to zero. In fact, for the present purposes there is no need to consider such a general 
limit and it is enough to consider that all the p^^, . . . ,p^^ that tend to infinity do so in the 
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same way, that is p^^ — p^^ — . . . — p^^^ — p. Therefore, in the numerator of the right hand 
side of Eq. (75) the greatest power of p [p^) will correspond to the numerator of the left 
hand side of the same equation. The same holds true for the denominators and therefore 
taking the limit p ^ oo will give the desired result. By using the generalized Gaudin's 
theorem we can write the trace of Eq. (75) in terms of the contractions 



pa 



(76) 



and therefore in order to evaluate the left hand side of Eq. (75) we will have to consider the 
matrix of contractions 



C 



pa 



lim C 

[p— >oo 



-^lim [(l + pT)-V] 



pa 



(77) 



To obtain the explicit expression for the above limit we will take into account the bipartite 
structure of the matrix T given in Eq. (35) as well as Eqs. (36) and (11) to write 



C 







-1 






)] 


1 












p j 



(78) 



The limit [p ^ oo] will be taken in two steps, first the probabilities p^ with u ^ /ii, . . . , /ik 
will be set to zero and the remaining Pn^jP^^, etc will be taken as equal to an unique 
parameter p that will be made to tend to infinity afterwards. By applying the first step 
(i.e. setting to zero the p^, {u ^ /ii, • • • , /^i^:)) the argument of the inverse matrix of Eq. 
(78) will be the sum of the identity matrix plus a matrix where the only surviving rows will 
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correspond to those of the probabilities p^^ going to infinity, that is 

/ ••• 



p 
p 



T-1 



X 



T vT-l 



X-' 



= I+p 







\ 







R/^kN T^j^i T^^2 







S^^N W'^il W^^2 



■ Sf,j^N W^Mifl ^Mk2 











/ 

(79) 

where we have set R = T = VX'^, S = Y'^X'^-^ and = X'^ to lighten a little 

bit the notation. Now it is convenient to introduce the vectors 



RfijN 



(80) 



and denote by e^. the Cartesian basis vector of dimension 2A'" along the direction ^j. The 
introduction of these quantities allow us to write Eq. (79) in a more tractable and compact 
form as 

K 

" (81) 



1 + pY, 

The expression of this matrix can be further simplified by introducing the 2N x 2K matrices 
E = (e^i • • • e^^e^^+N ■ ■ ■ (^iiK+N) and Q = {rfj,^ ■ ■ ■ Vij^^Sfj^^ • • • S/^k) in o^'^er to obtain the final 
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result 

1+ (82) 

The inverse matrix can be computed with the help of the Woodbury formula [10] 

{l + pEQ'^y^ ^I-pE{l + pQ'^E)'^Q'^ (83) 

where now the matrix to be inverted (I + pQ^E^ is a matrix of dimension 2K x 2K instead 
of 2A'" X 2N. With the p ^ oo limit in mind we will expand the right hand side of the 
previous equation in a power series of the inverse of p as 

{l + pEQ'^y^^I-E{Q'^E)-^Q'^ + ^E{Q^Ey^Q^ + --- (84) 

In a first sight one could think that this quantity, when multiplied by the remaining matrix 
depending upon p, will lead to a divergent quantity after the p — > oo limit has been taken. 
Fortunately, this is not the case and in order to recognize how different terms cancel out it 
is convenient to use the following identity 

+ pEE^ I ° ° I (85) 

that is only valid after the limit p,^ ^ {u ^ /ii, . . . , /ik) has been taken. The product 

{l + pEQ'^y^ [ ° I now simplifies owing to the fact that {l - E{Q'^ E)-^Q'^) EE^ = 

EE^ — EE^ — 0. Combining all the previous results together leads to the final result (with 
the limit p — > oo) already taken 






1 




+ E{Q^E)-'Q'' 



a (86) 



with = E^ \ j — i I . This result is quite relevant as it synthesizes in just 

ly yo 

one very compact formula the combinatorial number of contractions needed for the evalua- 
tion of the multiquasiparticle overlap of Eq. (86). The evaluation of the multiquasiparticle 
overlap could of course be carried out by means of the standard GWT but it would be 
a very painful procedure as the number of terms to be considered increases dramatically 
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(combinatorially) with the number of multiquasiparticle excitations (for instance, if we con- 
sider a four quasiparticle excitation overlap of a two body operator A (this is not such an 
uncommon overlap, see [8] for more details) then we should consider a twelve quasiparticle 
matrix element that involves 11!! (that is 10395) contractions instead of just the four need 
by using Eq. (86) . The above formula also presents computational advantages if one has 
to consider a variety of multiquasiparticle excitations. The reason is that it allows to split 
the calculation of the contractions in two well differentiated steps; the first is common to all 
the multiquasiparticle excitations to be considered and consist of the evaluation of the first 
matrix of Eq. (86). This is the most expensive computation from a computational point of 
view as it involves the inversion of the matrix X that is of dimension N x N. The second 
step depends on the indices of the multiquasiparticle excitations considered and consists 
of the construction of the matrices E and Q (of dimension 2N x 2K) out of the rows of 
the corresponding matrices and the inversion of the 2K x 2K matrix Q^E. For instance, 
in the case of a four quasiparticle excitation K — A and the later matrix to be inverted 
is of dimension 8x8. At this point, we can ask whether it is possible to generalize the 
multiquasiparticle overlap of Eq. (75) to the most general situation where the indices on the 
right multiquasiparticle excitation differ from the ones of the excitation acting on the left, 
namely 



(87) 



The answer is affirmative and involves the introduction of an operator T | ^ | that trans- 
forms the indices u into the indices /i in the following way 

■ ■ ■ I ^ I = (01 ■■■P.^. (88) 

By introducing this transformation operator we can express the overlap of Eq. (87) in the 
more familiar form of Eq. (73) 

~ ~ ~ -.-^ (01 • • • f ^ I 'T-'AP,, ■ ■ ■ hJ,. 10) 

PurPu2 ■ ■ ■ PvK^PtiK • • • P^i2P^ll 10/ \ 







I^Uil^U2 ' ' ' Pl^K^tiK ' ' ' Pn2Pni |0) 



■ ■ ■ I I ^-'P^K ■ ■ ■ 10) 



(89) 
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and therefore the right hand side of this expression can be evaluated by means of the 
contractions of Eq. (86) by substituting the matrix T representing T^^ by the matrix 

T I I representing T j ^ | T~^. The existence of T j ^ | , its explicit form as well 

V) VI VI 

as the expression of the matrix representing this operator are considered thoroughly in 

■ (A 

Appendix C. As can be observed in that appendix, the matrix representation of T I is 

VI 

nothing but a transposition matrix which exchanges given rows or columns of the matrices 
applied to it. Therefore, the matrix T j ^ \ representative of T j ^ | can be obtained 

VI VI 

very easily out of the one representative of T ^ by exchanging the appropriate rows or 
columns. 



VI. CONCLUSIONS 



By considering a certain limit of the statistical density operator where it becomes the one 
of a pure state we have been able to obtain the generalized Wick's theorem out of the corre- 
sponding statistical version (Gaudin's theorem). The advantage of deriving the Generalized 
Wick's Theorem (GWT) in this way is because the statistical version is much easier to derive 
and handle as it is based on the cyclic property of the trace over the Fock space. Using the 
limiting procedure we have been able to obtain the most general contractions needed in the 
GWT in a very easy way. By generalizing the limiting procedure we have also been able to 
obtain a GWT for multiquasiparticle overlaps. The corresponding contractions are given by 
simple and compact expressions that can accommodate easily many quasiparticle excitations 
reducing in this way the combinatorial complexity of the evaluation of multiquasiparticle 
overlaps to just a degree of complexity associated to the number of bodies of the operator 
whose multiquasiparticle overlap is required. The expressions obtained for the overlap are 
new and they will be very helpful in reducing the complexity of forthcoming beyond mean 
field calculation. 
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Appendix A: THE EXPONENTIAL OF ONE BODY OPERATORS AS CANON- 
ICAL TRANSFORMATION OPERATORS 

In this appendix we will remind the reader about the expression of exp(^)Q;pexp(— ^) 
where K — ^ ^nuO^nOiu is an one body operator written in terms of the matrix K and 
the fermionic operators of the condensed notation. Using the canonical anticommutation 
relations {a^, a^} = <J^u is very easy to verify that [K, ap] = — ^^((TKyi)pj,Q;,^ where the 
skew-symmetric matrix Ka = |(K — K^) has been introduced. Using the same rules we 
obtain ctp]] = X^,^(crKyi)^^Q;,^, etc which allows us to finally write 

eMK)apexp{-K) = + [it, a J + ^[k[k,ap]] + ... = J](e-'^^^)p,a, (Al) 

1/ 

The matrix M =e~'^^'^, owing to the skew-symmetric character of K4, satisfies M-^ = e^^"' ~ 
ae'^^^a or McrM^ = a that is nothing but the condition for the matrix M of being the matrix 
of a canonical transformation. 

Now we will show a result needed in the developments of the paper that states that if 
the matrices M and T both satisfy a relation of the type McrM-^ = a (i.e. they are matrices 
representing canonical transformations) then the relation 

(1 + MT)-V + (7(1 + T^M^)-^ = (7 (A2) 

holds. We start by considering 

a{l+T^Wf)-^ = ct(1-T'^M^+. . .) = a-T-^aWf+. . . = a-T-^M'^a+. . . = (l+T-^M-^)-V 

(A3) 

The right hand side can be written as (1 -|- MT)~^MT(j — a — {1 + MT)~^cr and from here 
the sought result of Eq. (A2) easily follows. 

Finally, by means of an example, we will argue that not all possible canonical transforma- 
tions satisfying McrM-^ = a are necessarily given in terms of a skew-symmetric matrix Ka by 
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the expression M —e To this end just consider the case of a bidimensional configuration 

/ a 

space. In this case, the most general skew-symmetric matrix is given by Kyi = I 

-K e" \ 

where a is a complex number and therefore e = . On the other hand, the 

y e-'' j 

most general bidimensional matrix M satisfying the canonical transformation condition is 



1 m ] 




( m\ 




1 and M = 




\ 1/m j 




[ 1/m / 



given by both M = and M = . Obviously, the first matrix is of 

\^ 1/m y ' ' 

the e~'^^^ kind but not the second one. 

Appendix B: CALCULATION OF THE DERIVATIVE OF CERTAIN TRACE 

In this appendix we will compute the following derivative 

-^IV[ln(l + MT(A))] (Bl) 

OA 

where both T(A) = exp(— AcrS) and M are matrices of respective canonical transformations 
(i.e. § is an skew-symmetric matrix). The logarithm in the previous expression has to be 
interpreted as the corresponding Taylor series expansion 

2 3 4- 

^jC ^jC 

ln{l + x)^x-- + --- + ... (B2) 
Now we have to consider the derivative of MT(A) with respect to A 

A (MT(A))^, = - (MT(A)a§)^, (B3) 

its square 

A {MT{X))l = A ^ ((MT(A))^ (MT(A)),,) (B4) 

T 

= - ^ ((MT(A)a§)^, (MT(A)),, + (MT(A))^, (MT(A)a§),,) 

r 

= - (MT(A)cr§MT(A) + MT(A)MT(A)a§)^^ 

and higher powers 

(MT(A))^^ = - (MT(A)(7§(MT(A))VmT(A)MT(A)(7§MT(A) + (MT(A))^MT(A)(7§) ^ 

(B5) 
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The generalization to higher powers is fairly simple and we only have to be careful with the 
non commutativity of the matrices involved. The results obtained so far for the derivatives 
are not very useful due to the increasing number of terms but fortunately we only need to 
consider its trace. Using the cyclic invariance property of the trace we can rearrange all 
different matrices in the derivatives to end up with a general and compact expression 

-^Tr [(MT(A))"] = -nTr [(MT(A))" a§] (B6) 

Using now this result we can write 

J-IY[ln(l + MT(A))] = ^Tr[(MT(A))]-l^iv[(MT(A))V^^^[(MT(A)f]^^^^^^ 

= -Tr [(MT(A) - (MT(A))^ + (MT(A))''^ + . . .) a§] 
= -Tr [(1 - (1 + MT(A))"^) dS] 



that leads to the final result 



-^Tr[ln(l + MT(A))] = -Tr [aS] + Tr [(1 + MT(A))"^ aS] (B8) 



Appendix C: THE OPERATOR OF INDEX TRANSFORMATION T 



/i 



The operator T | ^ | is defined is such a way that it transforms the quasiparticle an- 
nihilation and creation operators with index // into the ones with index u and viceversa, 
i.e. ^ 

/3^ = r |^^j/?,t|^^ I (ci) 

or 

/3^f l^^j =r |/?, (C2) 

for the quasiparticle annihilation operators (3/^ and (3^- It also has to leave the other quasi- 
particle operators unchanged, i.e. 

r\^\f3, (C3) 
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for any index a different from fj, and u. Finally we will impose a normalization that makes 



it to fulfill r 



operators P^. This normalization is not the natural one T 



where |0) is the vacuum of the quasiparticle annihilation 



\(f)) — — 10) (i.e. exchanging 



two fermionic quasiparticle states should lead to a minus sign) but is more convenient as it 
makes unnecessary to keep trace of the minus sign in the related expressions. However, if 
the reader feels more comfortable with the later normalization just multiplying by a minus 



sign the definition of 



/i 



below is enough. Using the notation of previous sections, one 



can write the requirement of Eq. (CI) and Eq. (C3) as well as the ones corresponding to 
the creation operators as 



/i 



/i 



r 



an 



pc 



where the bipartite matrix T | ^ | is block diagonal and given by 



T 



\x \ \ T{ii 1/) 

f{n ^ v) 



The matrix T(/x ^ i/) is a permutation matrix which has the structure 

/ 1 : \ 



T{fx ^ i/) = 



1 

: 1 : 



! 1 ! 

1 ••• 

! ! 1 



1 7 



(C4) 



(Co) 



(C6) 
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with all the matrix elements not explicitly given equal to zero. The matrix elements of 
T{ii <-> 1/) can be written in a very compact form as T{ii <-> 1/)-^ — 5ij + {Sii, — 5ifj){5ji^ — 5ji). 
This is a permutation matrix because, when applied to an arbitrary matrix A to the right, it 
exchanges its and v rows. On the other hand, if the matrix is applied to the left then the 
columns fi and v of A are exchanged. The permutation matrix T(/i u) is real, symmetric 

and idempotent (T'(/x <r^ v) = 1). It is possible to write the permutation operator T j ^ 

as the exponential of an one-body operator 



r 




exp(i^ 




(C7) 



with the one body operator X | ^ | = iX^/xi/^ ( ^ I OipOic- Using the general result 

pa 



V 

with T = exp(— crK) and taking into account that 



(C8) 




exp 




(C9) 



a little of algebra leads to the final result 



r 



V 



exp 



exp 



(CIO) 

Just for simplicity we will assume in this appendix that the creation operators are the 
hermitian conjugates of the corresponding annihilation ones and use the standard represen- 
tation for them. This assumption is not crucial and all the results presented in this 
appendix are independent of it. Now let us assume that we want to evaluate the matrix 
element . . ./Sj^^T^O/?^^ . . where Tq is an arbitrary operator carrying out a 

canonical transformation. In order to apply the main results of the paper we have first to 



transform {4'\f3yi^ ■ ■ ■ into • • • f^ix-^T \ ' " | by means of a transformation 

Vl ... Vk 

operator that in a first sight could be thought to be expressed as the product of elementary 
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transformation operators 

f - I =T| - ITI - I ...T| - I (Cll) 





However, a little care is needed as in the process of moving T | | to the left 



(C12) 

it might happen that one of the remaining quasiparticle operators have an index /ij that 

coincides with and therefore the transformation operator T I I turns it into /Xi and 

this will definitively alter the rest of the operators needed to carry out the transformation. 
An economical way to eliminate this possibility is to "discard" in the transformation operator 
those indices which are equal in both the fij and sets making it impossible that one of 
the elementary operators is changing more than one index. The mathematical expression 
for the operator would then be 

... Vk ) r=l,...K;ni^uj \^r) 

The order of the elementary operators is irrelevant as they commute among themselves 
and the only remaining detail is that the final product of quasiparticle operators (3^^^ . . . (3^^ 
contains the same quasiparticle indices as P^^^ ■ ■ ■ Pui but not necessarily in the same order. 
As the annihilation quasiparticle operators anticommute among themselves the reordering 
of the product to bring it to the desired form will introduce an additional sign that will be 
denoted fi, and is given by minus one to the number of transpositions needed to reorder the 
indices. This sign does not show up in the evaluation of the overlap matrix element 




{<l>\(3,^...(3,,TcP+...(3, 



Ml ■ ■ ■ ^Mjc I ^/ , , - / Ui . . . LtK , ^ 



Ul ... Vk 



(C14) 
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as it is the same both in the numerator and the denominator. It has to be kept in mind, 
however, that the sign has to be considered in the evaluation of the norm overlap (the 
denominator of the above expression). 



Finally, let us consider the combined effect of T 



(we now go back to a more 



compact notation for the indices /i and u) and an operator Tc carrying out an arbitrary 
canonical transformation with transformation matrix T {apTc — Yla ^pcr'^c^^a) 



apT 



fc = E^ 



OLaTc 



with 




(C15) 



Taking into account the decomposition of T 
mations of Eq. (C13) we can finally write 



T=i n 

^r=l,...K;ij,i^Uj 




(C16) 



as the product of elementary transfor- 




(C17) 



or in other words, the matrix T is obtained from the matrix T by the exchange of the 
rows Vr into the [ir ones. As an illustration of the procedure consider the transformation 
of l3^(3gl3i5(3-2(3Ql3il32i into /3i9/3i3/36/322/334/3i5/3i. The repeated indices are 6 and 15 and they 
will be discarded. Then the index 5 will turn into 19, the 9 into 13, the 2 into 22, the four 
into 34 and finally the 21 into 1. In this way we will end up with /9i9/?i3/?i5/?22/96/^34/^i and 
we will need four transpositions to bring it into the desired order (that is, fv — +1). The 
matrix T will be obtained out from the matrix T by exchanging its row number 5 with its 
row number 19; its row number 9 with its row number 13 and so one. The whole procedure 
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can be very easily implemented on a computer procedure. 
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